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ABSTRACT 


This  paper  presents  a  general  expression  for  the  WKB 

approximation  to  the  propagator  corresponding  to  an  arbitrary 

Hamiltonian  operator  H.  For  exa.npt  e,  if  the  correspondence  rule 

used  to  pass  from  the  classical  Hamiltonian  H  to  H  is  such  that 

c  ~ 

it  associates  aP.Q'I  +  (l-a)Q",p.  to  p.q-*,  then  the  formula  gives 
*•  r  *•  />  1  i 

\kB  "  Kwexpl(j5~a)/^  Hc/‘>qlapi)  (qc(t>  'Pc(fc)  'fc)dt}  * 

/y  "T  L 

where  =  (2itifi)  /  (det  M)  ^  exp(iSc/h)  is  Van  Vleck's  well- 

known  formula,  S  being  the  action  functional  evaluated  at  the 
c  2  i  i 

classical  path  (q^.p^)  and  sc^<ia  More  generally, 

the  formula  presented  here  applies  to  any  system  described  by 
a  function  f(x,t)  whose  time  evolution  is  given  by 
(H(x,  k'i/dx,  t)  +  k^/it)f(x,t)  =  0,  regardless  of  the  form  of  H. 
The  Schrddinger  equation  of  quantum  mechanics  and  the  Fokker- 
Planck  equation  of  diffusion  are  obvious  examples.  Many  examples 
are  discussed.  This  generalizes  results  obtained  in  a  previous 
publication  (J.  Math.  Phys .  1_8  (1977),  786-90). 
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INTRODUCTION 


The  purpose  of  this  paper  is  to  obtain  a  general  expression 


for  the  WKB  approximation  to  the  propagator  corresponding  to  an 

.  1 - - 

arbitrary  Hamiltonian  operator.  In  an  earlier  publication  we 


<T 


determined  a  range  of  validity  of  Van  Vleck's  well-known  formula  , 

3 

which  was  known  not  to  be  universally  valid  .  The  approximation 
derived  here,  generalizing  Van  Vleck's  formula,  is  also  valid  beyond 
quantum-mechanical  applications^as  it  applies  to  any  system 
described  by  a  function  f(x,t)  whose  time  evolution  is  dictated  by 

£  H(x,  k*>/dx,  t)  +  kVStJ  f(x,t)  =  0/  (1) 

such  as  the  Fokker-Planck  equation  for  diffusion  processes. 


THE  GENERAL  WKB  APPROXIMATION  FORMULAE 


^We  operate  in  n  dimensions  and  summation  over  repeated 
indices  is  implied.  The  following  theorem  summarizes  our  findings: 


Theorem 


Let 


•  H (Q, P, t)  be  an  arbitrary  Hamiltonian  operator  (the  lack  of 

constraints  imply  that  it  could  be  non-Hermitian,  time- 
dependent,  non-quadratic  in  P,  etc.) 

M 

•  H  (q,p,t)  =  H(Q-*q,  P-»p,  fi  =  0)  be  its  classical  counterpart. 

c  «>  ~ 

•  The  correspondence  rule  used  to  pass  from  H  to  H  be  such  that 

C  A, 

it  makes  the  following  associations:  ’ 

f(q)  f(Q),  f  {?)«-»  f(P)  (2) 

p.q3<L«>  ap.Q-5  +  (1  -  aX^P.  for  a  given  a.  (3) 

(Note  that  (2)  and  (3)  do  not  imply  Hermiticity  of  H, 
even  if  a^h) . 


•  (q  (t),p  (t)  )  be  the  classical  solution,  solving  Hamilton's 
c  c 

equations  for  H  such  that  q  (t  )  =  q  and  q  (t,  )  *»  a  , 

r»  a  *  >  n  ' 


•  VVWV  5/fPci(t)C(t)  "  Hc(qc(t)'Pc(t,'t)]  dt  (4) 

be  the  classical  action  (  Ts  [  T  £ 

.  M1.(qb,tb.qa.ta)  =  -VV^X3  (5) 

be  the  Van  Vleck  -  Morette  matrix,  with  determinant  det  Mj 

•  K(q,  ,  t,  ,q  ,t  )  be  the  propagator  corresponding  to  H,  defined  by 

b  b  a  a 

f H{Q,P,t  )  -  ifi  1 K  =  0,  (6) 

lim  K  -  $  (q, -q^)  ,  (7) 

t  -¥  t  »  a 

b  a 

where  Q  is  represented  by  q^  and  P  by  -ih 

•  (q^,  t^,  qa>  ta)  be  the  WKB  approximation  to  the  propagator, 
defined  by: 


1Cb[“(S’.p'V  -  if^‘b]  Vb  - 


Z  (q^-q _>. 

b  a 


Then  the  WKB  approximation  is  given  by 


vb 

"Wb  ’  %V  exp[  (!i  -  a)  f  lil-C  (qo(t)'  Po<t)'  *)  dtl  .  (10) 


where  K  is  Van  Vleck's  formula: 
W 


■  (2irifi)  n (det  M)*5  exp(iSc/fi). 


(The  case  where  det  M  =  0  is  not  examined  here) . 
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More  generally,  the  WKB  approximation  is  given  by 


WKB 


A  exp(iS  /fi) 
o  c  ' 


where 


»A  -  K  (0,  t  ,  0,  t  )  , 
o  o  b  a 


(12) 


(13) 


•  X  (q,  ' ,  t,  ,q  '  ,t  )  is  the  propagator  corresponding  to  the 
o  o  oa  a 

Hamiltonian  ooerator 


H  5  Jjglj  (t)u  +  Hf  (t)vlj  +  k. 3  (t)w1 .  , 
o  ~  13  1  ~  J 


(14) 


^pPfj \yv» 


f . . (t)£ 

/  13 


viw 

(15) 


re*1*) 

The  correspondence  rule  used  to  pass  from  to  H  is  such 


that  it  makes  the  following  associations: 

PiPj*-»Uij  q1q^  V1^  q^p,  W1 

■*»  /V  '  J  <V  J 

PROOF  OF  THE  GENERAL  WKB  APPROXIMATION  FORMULAE 


(It) 


Let  us  begin  by  giving  a  simple  example  illustrating  the  fact 
that  the  Van  Vleck  formula  in  (11)  is  not  always  equal  to  the 

WKB  approximation  Let  H  be  an  operator  such  that  -  K^. 

The  operator 

H'  5  F-1(Q)HF(Q)  *  (17)  _ 

has  the  same  Hc  (and  hence  the  same  K^)  as  H.  Yet  its  WKB  appro¬ 


ximation  is  not  K  but 

W 


-1 


Kkb  *  F(q.,F  tqb>V  - 

as  can  be  verified  by  direct  substitution. 


(17a) 


f 


We  now  prove  the  theorem.  Formula  (10)  will  be  proved  by 
generalizing  the  proof  given  in  ref.  1.  The  more  general  formula 
(12)  will  be  proved  by  generalizing  the  method  used  in  ref.  4, 
which  involves  path  integrals.  It  is  presented  separately  because 
it  is  more  on  the  heuristic  side. 

Proof  of  first  formula 

In  ref.  1,  where  we  investigated  the  range  of  validity 
of  Van  Vleck’s  formula,  we  assumed  that  the  WKB  approximation  was 


of  the  form: 


\kb  =  a^c  exp(iB/fi)/ 


where  C  and  B  are  real  functions  of  q„  ,  t,  ,  q  ,  and  t  ,  inde- 

o  d  a  a 

pendent  of  fi ,  and  found  that 

..  *(  /  .i  ..  Cl  \  i/  .  . 


-i  +  rmn.(B',,ht,  )+»[.(£ 

'  (sfT^S  '•>)„,,('”<'»* 1 TO>) 

-^’iaS-TOlf-lTOW,^ 

1  missis] 


01) 


where  B"  5  >2B/fcq,  1 'fcq.  ^  ,  etc.,  and  F  is  Cohen's  F  function. 
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establishing  the  correspondence  between  H  and  H  by 

c 


H 


dp  dg  du  dv  F(u,v)H  (q,p,t) 

c 


(2r*)~2n  fc 

t+  J  w  ^ 

exp  £  (i/fi)£(q-Q)  .u  +  (p-P).v3j  *  (20) 

P  effects  the  generating-function  correspondence 

F(u,v)  exp  £  (-i/fi)  (Qu+Pv)J  exp£(-i/ft)  (qu+pv)3  ,  (21) 

and  a  set  of  Fs  effecting  a  given  correspondence  H  4+  H  can  be  found 
^  c  ~ 

by  solving 

F(u,v)  fdp  dq  H  (q,p,t)el(qJ+pv)/Jfl  =  (2rfi)  tr  (e1  (2U+~V)  ^H)  .  (22) 


r(Rtn 

This  equation  for  F,  to  be  understood  in  the  sense  of  distribution 
theory,  may  or  may  not  have  a  solution,  and  the  solution  may  or  may 
not  be  unique  if  it  exists. 

In  ref.  1,  F  was  assumed  to  be  a  function  of  u.v/h  only, 
and  this  requirement  is  maintained  here,  as  u.v/fi  is  the  only  dimen¬ 
sionless  quantity  one  can  form  with  u,  v,  and  and  it  is  assumed 
here  that  H  contains  no  constants  (other  than  -fi)  that  do  not  appear 

~  7 

in  Hc«  It  can  be  shown  that  the  fact  that  F  depends  on  the  combi¬ 
nation  u.v/fi  implies  that  the  operator  H  corresponding  to  pmqn  in 
1  dimension  is  a  linear  combination  of  all  the  possible  ordered 
arrangements  of  P  m  times  and  Q  n  times  (true  "factor  ordering"). 

**  ~  g 

Thus,  no  true  divisors  are  allowed  in  Cohen's  scheme  ,  an  important 
restriction  since  the  Laplacian  in  curved  spaces  contains  true 
divisors.  We  shall  return  to  this  point  later. 

9 

In  ref.  1,  F  was  also  assumed  to  be  real  so  as  to  insure 
the  Hermit icity  of  H.  This  requirement  is  dropped  here.  Therefore, 
equating  both  the  real  and  the  imaginary  part  of  the  constant  term 
and  the  term  proportional  to  fi  to  0  in  (19)  yields  the  following 
4  equations  to  be  satisfied  by  B  and  C: 


Re) 


Hc<  V 

2J? 

H  (  v 

'*&  , 

c  b 

n 

1  Me 

2  'ii 

L 

.  tb)  *e[F(0)]+  |S  =  0 


(23) 


(24) 


P  53  *>  S/j>q, 


J 


3L  (C2)  +  2C2  (q,p,t  )1 


im 


f^r 


(25) 

[F'  (0)]  =  0 


p  =^B/aq, 

o 

q  "  <Iv 


jRe[F'(0)]-  H  Im[F(0)3^  f*  .  >* 

L  C  ‘P  = 


'*B/Zq. 


(26) 


q=q^ 


Im[F(0)]  T  c['  ( q^ »  p <  tb)  +  H  3V.  ^  1 

L  C  ^  p^3/>% 


o. 


where  the  third  equation  was  rewritten  in  a  more  compact  form. 

First,  let  us  translate  the  conditions  (2)  and  (3)  on  H 

A# 

into  conditions  on  F.  These  are: 

F  (0)  =  1;  Re[F'(0)]=  0,  Im[F’(0)]=  a  -  h  .  (27) 

Indeed,  (2)  readily  results  from  (21)  with  F(0)  =  1.  If  we  diffe¬ 
rentiate  (21)  with  respect  to  v,  then  set  v=0,  then  do  the  same  with 
u,  we  obtain 

pq  <r>  F(0)  (QP+PQ)/2  -'fF’(O),  (23) 

which  yields  the  2  conditions  on  F’(0). 

It  is,  in  fact,  possible  to  drop  the  requirement  that  the 
correspondence  be  given  by  an  F  function,  so  that  the  results  are 
valid  for  any  H  satisfying  (2)  and  (3).  This  will  be  seen  when  we 
treat  the  more  general  formula  (12). 
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Let  us  now  solve  (23)- (26)  for  B  and  C.  Since  F(0)  =  1, 

(23)  and  (24)  are  recognized  as  being  one  of  the  Hamilton- Jacobi 

equations,  yielding  B  =  S  .  Note  that  since  (23)  and  (24)  stem 

c 

from  equating  only  the  constant  term  to  0  in  (19),  an  approximation 

to  zeroth  order  in  fi  yields  B  -  S  and  no  information  on  C  other 

c 

than  the  fact  that  it  must  be  such  that  the  boundary  condition  (9) 
is  satisfied  (this  is  not  sufficient  to  fix  C)  . 

To  solve  (25),  we  note  that  if  Xm[F'(0)}=  0  (a=h)  then 

2 

C  =  det  M,  since  the  (continuity)  equation  satisfied  by  the  Van- 
Vleck  -  Morette  determinant  is  precisely 


3.  r 


]}c  <  V 


det  M 


1  +  3L  Pdet  M 1  =  0  . 


This  leads  us  to  write 


(29) 


C 


2 


=  N  det  M, 


(30) 


which,  when  substituted  in  (25)  and  using  Im^F' (0)}  »  a-h  yields 


the  equation  for  N(q,  ,t  ,  q  ,t  ); 

D  D  cl  S 


■*v 

with  N*^l  as  t,->  t  and  q,- 

OS  3 


+  N(2a-1) 


-i) 


■%  if  KW 


83  0 . 

5(qb-qa)  as  tb^ta 


(31) 


The  search  for  a  solution  of  (31)  is  facilitated  by  the 
observation  that  the  first  two  terms  form  the  convective  derivative 
of  N  with  respect  to  the  final  endpoint.  Knowing  that  the  convec¬ 
tive  derivative  of  any  function  f  of  position,  momentum,  and  time 
(at  the  classical  path)  with  respect  to  either  the  initial  or  the 
final  endpoint  is  0,  i.e., 

LV 

(to  be  proved  shortly),  we  are  led  to  a  trial  solution  of  the  form 


q  3 


c  b  n 


-1 

"V-1 


f (q  (t),p  =  o 

C  C 


(32) 
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N  *  exp 


£  cx-a.,  f  2^1 

1  x  ^.v 


(q  (t),p  (t)  ,  t)  dtl  . 

c  c  -J 


By  direct  substitution  in  (31),  (33)  is  seen  to  be  the  correct 

solution  (use  f=q  and  f=p  in  (32)). 

c  c 

To  prove  (32),  we  first  observe  that  it  is  sufficient  to 

prove  it  for  f=q  and  f=p  :  the  chain  rule  will  then  extend  its 
c  c 

validity  to  any  f.  Now,  derivatives  of  the  classical  solution  with 
respect  to  any  parameter  introduced  by  the  boundary  conditions 
(here,  t  ,  t  ,q  ,  or  q  )  are  known  to  be  solutions  of  the  equation 
of  small  disturbances,  obtained  from  the  second  variation  of  the 
action  functional10  13 .  Since  this  equation  is  linear,  its  general 
solution  is  a  linear  combination  of  2n  linearly  independent  solu¬ 
tions.  Thus,  the  only  solution  vanishing  at  both  t  and  t  must 

3  D  • 

be  zero  everywhere.  Now,  the  left-hand-side  of  (32)  with  f=q^  (t) 

is,  by  its  very  construction,  a  solution  of  the  small-disturbance 

equation.  It  vanishes  at  t  because  q  (t  )  =  q  (a  constant).  It 

a  c  a  a 

also  vanishes  at  t  because 

rsc'l  ,-;;<«*£  (W  .... 

(since  7)  /  ^  t^  commutes  with  <V*t  when  acting  on  q^ (t, t^, q^) ) . 

Thus,  (32)  is  true  for  f=q  .  The  case  f=p  is  proved  by  observing, 

c  c 

by  substitution  in  the  equation  of  small  disturbances  in  phase  space 
(O'  h  =  0)  ,  that  if  u  is  a  parameter  introduced  by  the  boundary 

A* 

conditions,  then 


'K  ~  C  1)*. 
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where 


k+d/dt\ 


[  ^Pc(t)Au  ) 


(34) 


-1  A 

D  s  g  ( 1  f.  -  k), 

and  f,g,  and  k  are  defined  in  (15).  This  completes  the  proof  of 

the  WKB  approximation  formula  (10)  . 

Note  that  (19)  indicates  that  when  H  is  quadratic  in 

*.  c 

both  p  and  q,  the  term  proportional  to  fi  is  0  because  C  is  inde¬ 
pendent  of  q  and  g,  (S  being  quadratic  in  q  and  q,  ),  and 

3.  O  C  HO 

higher-order  terms  are  0  because  they  involve  third  and  higher 
derivatives  of  H  .  Thus,  the  WKB  approximation  is  exact  in  that  case. 
This  goes  beyond  the  well-known  result  because  H  does  not  have 
to  be  Hermitian,  so  that  the  extra  exponential  term  in  (10) 
supplementing  the  Van  Vleck  formula  is  not  constant. 


proof  of  second  formula 

Formula  (12)  will  be  proved  using  the  path- integral 
approach,  generalizing  a  method  presented  in  Ref.  4.  More  details 
on  this  method  and  its  extension  to  a  WKB  expansion  of  the  pro¬ 
pagator  for  arbitrary  Hamiltonians  will  be  presented  elsewhere. 

The  propagator  K  can  be  written  as  a  phase  space  path 
integral  as  follows: 


K<VVW 


■  Ij^)  exp  ( iSAO, 


(37) 


where  S  =  f [pq  -  H(q,p,t)J  dt  is  the  action  functional  and 

J*  is  the  space  of  paths  (q,p)  such  that  q(tQ)  *  qa  and 

q(t  )  “  q,  .  If  S  is  expanded  around  the  classical  path  (q^.p  ), 
b  *3  ”  c 

its  first  functional  derivative  vanishes  by  definition  of  the 

classical  path  and  we  obtain 


K 


iSc/ft  C  r  <^x^y1p,*-S^(x,y) e  i-0-c(x,y)/ft 

L  l.”  J  ' 


(38) 


where  S  is  the  space  of  paths  (x,y)  such  that  x(t  )  * 
o  _  a 

x ( t  )  =0,  J \L  contains  the  terms  beyond  the  second  functional 
b  c 


derivative,  and  the  second  functional  derivative  is: 

S£(x,y)  *  y'dt^yi(t)x1(t)  -  bg1-*  (t)  y  i  (t)  y..  (t)  -  Jjf  _  (t)x1  (t)x-*  (t) 


~  k1j  (t)y.(t)xj(t)] 


with  f,  g,  and  k  in  (15). 


We  can  define  a  measure  w^  normalized  to  1  and  absorbing 
the  second  variation  of  S  by: 


dw  (x,  y)  5 


f  tA  x  <A' 


■  t 


[J  exp(iS^(x,  y)/b) , 


the  normalization  factor  being 


A  - 
o 


f  [ilirl 

Jf.  L  h"  J 


exp ( is " (x, y)/h). 
c 


Now,  it  is  observed  that  the  S"  term  in  (39)  is  in  the  form  of 

c 

an  action  functional  corresponding  to  the  fictitious  Hamiltonian 
Hq (x,  y  »  hrjX^  (tjy^j  +  %£ (t)xXx^  -  k\  (t) y^  .  (42) 

Hence,  (37)  indicates  that  A  must  be  the  propagator  K  (q  ’,t  ,q  ’  t  ) 

O  O  O  D  a  cl 

corresponding  to  H  ,  evaluated  at  q  1  =  q  ’  =  0.  But  for  which 

o  3  /  d 

Hamiltonian  operator  H  ?  It  makes  sense**  that  it  should  be  the 

«-o 

operator  derived  from  Hq  using  the  same  correspondence  rule  linking 
H  and  H,  i.e.,  (14)  with  (16).  This  leaves  us  with 

C  /V. 


K  =  A  e  c 
o 


iS_/h  C  e-iil_(x,y)/h 

JsO  6 


dw(x,  y)  . 


It  can  be  shown  in  this  general  case,  as  was  done  before  for 
^  ^  13 

special  cases  ,  that  the  expansion  of  the  il  term  followed 

c 

by  the  evaluation  of  the  path  integrals  (the  correspondence  rule 
being  taken  into  account)  yields  a  series  in  f», 

K  =  A  eiSc//fl  (1  +  RK.  +  R2K0  +  .  .  . ) ,  (44) 

O  12 

which  identifies  the  constant  term  as  the  WKB  approximation.  ■ 
Let  us  now  retrieve  formula  (10)  from  this  more  general 


case.  The  operator  Hq  in  this  case  is 


12 


as 


HQ  2  **gXJ(t)P.P  +  Hf  (t)Q1QJ  +  k1  (t)faPiQ;j  +  (l-a)Q:}PiT/  (45) 

which  can  be  rewritten,  usinj  [?  'pil-  111  fV 

Ho  a  Hoo  +  ift(>s-a)k1i(t) ,  (46) 

where  ~  ^ 

Hoo  2  (fc)PiPi  +  15fi-i  (fc)Q1Q:’  +  l5k1.  (t)  (P .Q:’-K);iP.) .  (47) 

-  ^  'v  .  J  J  v  J  a.  ^  X 

v  v  *v  ^ 

Now,  since  H  is  quadratic  and  Hermitian,  its  propagator  K  is 
*<oo  oo 

given  exactly  by  Van  Vleck's  formula, (11) .  In  this  case,  however, 

the  S  in  (11)  is  zero  because  q  '  =  q,  '  =0.  (In  fact,  in 
c  a  b 

general  S  *  q,  1  Xp  .  (t  )  -  q  '  1p  .  (t  ) ,  where  (q  ,p  )  is 
c  b  rco ,  i  b  rco,i  a  Mco  co 

the  classical  solution  for  Hq) .  Further,  the  "det  M"  in  (11)  is 

the  same  as  the  "dct  M"  for  H  because  H  and  H  share  the  same 

c  o  c 

equation  of  small  disturbances  and  M. .  is  a  boundary  value  of 

ii-13 

a  specific  solution  of  that  equation  .  Therefore,  K  ■ 

-n/2  h  00 

(2rrifi)  (det  M)  .  Now,  if  h  =  H  +  f(t),  the  propagators  K. 

~1  ~  2  _  1 

and  K  are  related  by  K,  “  KL  exp(-  (i/fi)  J  f(s)ds).  H  and  H  are 
2  12  x  ^o  «oo 

related  in  this  manner.  Putting  all  these  results  together  gives: 


-  (2ffift)  n/2(det  M)**  exp  |  i 


S  + 
c 


(*S-a)  J  k1i(s)ds  | 


(48) 


which  is  formula  (10) .  Note  that  this  suggests  that  the  norma¬ 
lization  factor  (2«i)  )n^  is  universal  and  independent  of  H 

c 

Note  in  passing  that  it  is  not  always  easy  to  find  out 

what  p.qJ  corresponds  to  given  the  correspondence.  Scaling 

X  c  ** 

tricks  (replacing  Q  by  A Q  in  functions  of  Q,  then  differentiating 
with  respect  to  A  and  setting  A  equal  to  0)  sometimes  help. 


SOME  EXAMPLES 

We  begin  with  an  example  pointing  out  that  formula  (10) 
is  not  restricted  to  the  correspondence  rule  being  effected  by 
an  F-function.  Consider  the  Hamiltonian 


H  =  *59  ^(Qjfp^  -  A  (Q)1  glS(Q)g1J  (Q)fp.  -  A  (Q)1g-!t(Q)  +  V (Q)  (49) 

corresponding  to  Hc  =  (q)[  p^-Ai  (qj][Pj  -  A..  (q)J  +  V(q),  where  (50^ 

g  g^  =  J. 3  and  g  -  det(g.  . ).  Since  it  is  Hermitian,  a****.  There 
lk  ~  x  ~  id  15 

is  no  F  in  general  because  of  the  divisors  .  Nevertheless,  a  direct 
substitution  shows  that  Van  Vleck's  formula  applies,  and  we  get: 

+  rA  **  f  ^ 

b  *  M  «  M  Zm2- 


-  H,;|  I 


r*  - r* r"  +  r*  .  )  (a) 
=  o(tf),  ( 


(s-i) 


where  the  following  properties  and  definitions  were  used: 


M=  det  M 


Qi  S  %ia<5ja.k  +gka.j  ~  3jk,a>  ) 


_ij  .  p4'  mj  pi  r,im  r*  (rz) 

9  ,k  ■  ■  U3  '  r„V  ,  9  ,i  ~2«1  ru>  ' 

9c3<V  -  ^X’lAct'V  -  WJ  ;  pci<V  m*sc/*%\ 

7)tb 

as  well  a3  (29)^. 

On  the  other  hand,  (19)  gives  the  following  expression  for 


the  miss  term; 


K"1^  -  ifi3/J>tjK  =»  gij  ..(?«F(0)  -  iF'(O)  -  F"  (0) ) 

w  y  b  w  L  ' ij 

p  * 

+gij  .  (F  (0)  -  2iF'tO))£i  +  gl3F(0)^JJ  *1  (53) 

•3  c  c  J 

an  expression  which  cannot  be  matched  with  (51)  for  any  F,  for 

H 

C  =  M  . 
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Consider  now  the  Fokker-Planck  equation  of  diffusion  processes 


r  i  -2 -  rDlj(q,t)p]  -  JL  [vX(q,t)p]  ,  (54) 

"Dtr  V  17 

where  D  is  the  diffusion  matrix  and  v  the  drift  vector  •  It 

formally  corresponds  to  a  Hamiltonian  H  =.  P. P .D1"5  (Q)/2ifi  +  P-vX(Q) 

~i~j 

with  classical  Hamiltonian  H  «  p.P.D  /2ifi  +  p.v  .  since  the  P 

c  1  j  1 

factors  precede  the  Q  factors,  a  is  simply  equal  to  1.  Thus,  formula 
(10)  gives; 


Kwrb  =  (2irtfO"n/2(det  M^expJ"^ 

-  hf v1  (q_i(t))dtl. 

X  ' 1  u  J 


S  -  (2ifi) 

K  c 


•v 


pci(t,D 


ID 


j (qc(t))dt 
(55) 


Note  that  the  dynamical  equation  gives  p  . (t)  =  ifiq  3 (t) (D  ) . . (q  (t)) 

Cl  C  J  X  c 

so  that  in  one  dimension  part  of  the  integral  can  be  performed, 
yielding 


Vs  -  KwtD(qb,/D('5a,r’i  e*P(,*jr«'®'-v-D)/DMf). 


(56) 


In  the  ca.se  of  the  backwards  equation  (Q  precedes  P,  a*0) ,  the 
factors  \  are  replaced  by  in  (56)  and  the  broket  of  (Vs”). 

For  constant  diffusion  parameter  (D^l)  and  linear  drift 
v  =*  -Vq,  one  retrieves  the  well-known  propagator 


18 


W  Kt 


=  _*  —  W- 


I  - 


]■ 


1*1) 


which,  H  being  quadratic,  is  also  exact  (it  satisfies  (54)  exactly). 

It  can  be  shown,  by  direct  calculations,  that  the  "miss 
factor"  for  (55)  is  exactly  as  given  by  (53)  with  F(x)  ■  exp(ix/2). 


15 


Another  interesting  application  is  the  "lognormal"  process 
with  Hamiltonian 

2  2 


H  =»  d  Q  P  + 


ft  QP 

-v  V 


(58) 
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in  one  dimension,  useful  in  modeling  population  growth  .  Using 

2  2 

formula  (10)  on  H^  »  c*  q  p  +  ft  qp  one  gets: 

Vb  =  (Vqa,1’**Kw  e*0[/5<a-'s'TJ  ,  <59> 

20 

where 


Kw  =  (2tjifi)  Js(2^Tqaqb)'JSe:. 


“p[*b(uJi-'*T)l3-W0) 

21 

In  our  case,  a=0.  Note  that  for  a=0  the  exact  propagator  is 

K  »  KWKBexp(-i-ft«iT/4)>  (61) 

so  that  the  expansion  of  the  exponential  gives  the  terms  of  a 
WKB  expansion  of  K,  useful  for  checking  general  formulas. 


Let  us  also  mention  the  elements  of  the  WKB  approximation  for 
Hamiltonians  in  one  dimension  of  the  form 


,  m  n 
Hc  =  kp  q 


(Cz) 


where  k  is  a  constant.  The  classical  equation  of  motion  is 

2  J2.  t!  _i 

q  ®  (n/m)q  /q  and  the  Lagrangian  is  L  =  (q/q)^'*  (mk)”*'  (1-m  ). 

c  c  c 


For  m  ^  n  these  elements  are; 


q  (t)  **  A (t-t  )** 

^c  o  / 


I 


P  (t)  =  [k(m-n)]  (t-t  ) 


YT 


-t 
w  I 

L 


M  -  ^Tk(m-n)]M“ 


t 

1  ° 
w\ 

h-yt-b 
\~Y  ' 

w 

(m-n),*wi 

(qb"’ 

-  q 

a 

1 

! 

'Vb1' 

-n/m 

(qb 

—VI 

—  q 

a 

**(  w\- \\ 

J3L.  t-o 
o'  / 

'  -  W 


L 


(63) 
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For  m  «  n,  these  elements  are:  ^ 

qc(t)  =  exp(At+B)  f  qc(t)pc(t)  =»  (AAm)" 

A  s  T^Log  (qh/qA)  ,  B  5  Log  [  q*  (^/%)  fc«/TJ 

Sc  *»  (mkT)  (1-m  }  j^Log  (q^/q^j 

zL  _|  2*  w> 

M-  (mkT)**  [  (m-Dq^a^fLog  (q^/qb)J  ^ 

In  all  cases, 

,  ,  n(Jj-a) 

KWKB  *  Kw  qb^a 


(IH) 

(£  5) 


CONCLUSION 

We  have  produced  a  formula  to  approximate  the  propagator 

corresponding  to  any  system  described  by  a  function  whose  time 

of 

evolution  is  given  by  a  partial  differential  equation  that  isAfirst 

order  in  time.  This  approximation  can  be  supplemented  by  correction 

terms  that  can  be  generated  by  path  integrals,  and  this  will  be  the 

22 

subject  of  a  follow-up  study 
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is  fulfilled  by  an  F  function  iff  (q)  +  *|(q)  *»  k  ^(q),  where 
k  is  a  constant,  (q)  =  ^sg1^  an<3 

"Id) (*- i )  r*  q  +  (2^  +±) Q  rj 

+i  r-f  r4.  -  J-  r*  - +x\  rM.  -i . 

Z  ,'4  ij  *f  -j;  «/  '"‘,1 J 
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example  for  k=*0.  In  particular,  there  is  no  F  in  general  for 
the  case  =0,  In  one  dimension,  ^  k^  for  all  g1^  °£  fc^e 
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C(Hb  -  ifi  2/ttb,Kw  ='fi2  +  °<vta>  +  o(V’a,3  ' 

where  r  h  g'j/ri- H  +  c 1 .  r? 

<1  C  *),'  vj,f  '  *"j  '  -  I  *>  «*•■*  l 

is  the  Ricci  scalar curvature  tensor.  In  this  case,  -a  small-time 

approximation  is  the  same  as  a  WKB  approximation. 
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^  2  ^  ^  a  A  ^  / 
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21.  Exact  propagators  for  Hamiltonians  of  the  form  H  **  Z-a.Q  P  in 

one  dimension  can  sometimes  be  found  using  the  identity 

[N  ,  •  _  w\  fj  , 

1  -  [la.  Xn 

J  I  ,■!  1  (n.;)!  J  '  ,, 

i  l 

where  9(x)  ■  0  for  x<0  and  1  otherwise.  Thus,  for  T  =  X  a-x  (d/dx) 

\s%  1 

one  obtains 

T  T  *  P^/  \  T  7  «  irt. 

e  £<lt>  =  «: f  ("'  ***  1  A,  ;  ”T7S:  J  • 


For  ^=2,  use  of  the  integral  representation 

2  a  2 

sn  , .  -*S  /  -x  /4s  +  nx  . 

e  =>  (4TTs)  I  e  dx  gives 

fR  2 

eTf(x)  **  (If tra  >  ^  f  dy  e  y  '/4a2  f  (xey+ai’1’a2)  .  Changing  variables  to 
2 

u  »  x  exp(y+a  -a  )  gives 

L  » 

eT  &  (x-x* )  =»  (x’)  1(4ma2)  **  exp(  -(Log(x’/x)  -w^-a^  2/4a2) .  Since, 
it  will  be  recalled,  the  propagator  K  is  equal  to 

^-iTH/ft  J(  q  )  wiien  h  is  time -independent,  the  method  can  be  applied 


here. 
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